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We perform an analytical study of the scattering matrix and bound states in problems with many 
physical coupled channels. We establish the relationship of the couplings of the states to the different 
channels, obtained from the residues of the scattering matrix at the poles, with the wave functions 
for the different channels. The couplings basically reflect the value of the wave functions around the 
origin in coordinate space. In the concrete case of the X(3872) resonance, understood as a bound 
state of D°D*° and D + D*~ (and c.cQ), with the D°D*° loosely bound, we find that the couplings 
to the two channels are essentially equal leading to a state of good isospin 1 = character. This 
is in spite of having a probability for finding the D°D*° state much larger than for D + D*~ since 
the loosely bound channel extends further in space. The analytical results, obtained with exact 
solutions of the Schrodinger equation for the wave functions, can be useful in general to interpret 
results found numerically in the study of problems with unitary coupled channels methods. 
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I. INTRODUCTION 

The A (3872) resonance, observed by Belle [l| and con- 
firmed by CDFII, DO and BaBar collaborations [H-IH, has 
been the object of intense debate from the theoretical 
point of view (see recent workshop on charm exotics at 
Badhonef (H). Although different tentative explanations 
to its nature have been provided [(l-fl2j the idea most sup- 
ported recently is that it corresponds to a loosely bound 
state of DD* I12l UJ-[27| or slightly unbound, virtual 
DD* state fl3L l28l| . However, the energy of the resonance 
is very close to the D°D*° threshold, with the eventual 
charged components D + D*~ bound by about 8 MeV. 
The binding of the D°D*° could be so small as to render 
the relatively very bound charged components irrelevant, 
at least from the probability point of view, given the fact 
that the loosely bound component would extend much 
further in space than the c harg ed components. This is 
the idea behind many works [TJ [H, HH-fH, . However 
it was found in [l3[ that the couplings of the resonance 
to the charged and neutral components were practically 
identical, implying a near I — nature of the resonance 
as experimentally established. A pure D°D*° component 
would have an equal admixture of / = and 1 = 1 and, 
according to [13] would produce a ratio of 



B(X -> J/Vtt+tt-) 



(1) 



much larger than experiment. Indeed, in [29j several 
works based on pion exchange as the source for the DD* 



1 Fiom now on, when we refer to D°D*°, or DD* we 

are actually referring to the combination of these states with their 
complex conjugate in order to form a state with positive C-parity. 



binding are analysed thoroughly, stressing the impor- 
tance of taking into account the neutral and charged 
components to properly study isospin violation in the 
A (3872). There seems to be a contradiction between the 
intuitive idea of a dominance of the loosely bound compo- 
nent and the fact that experiment demands clearly an im- 
portant contribution from the charged components. The 
clarification of this puzzle and establishing the meaning 
of the couplings in terms of wave functions is the purpose 
of this present work. While in the theoretical calculations 
one normally uses field theoretical methods to evaluate 
observables [12|], without resorting to wave functions, the 
clarification of the puzzle forces one to face this problem 
solving the Schrodinger equation for the wave functions 
in coupled channels. 

For simplicity, we will assume that the A (3872) mass 
is below both the D+D*~ and D°D*° thresholds. The 
work proceeds as follows: in the next section we make a 
brief summary of (l3j to expose the problem. In sections 
III and IV we solve the Schrodinger equation in the case 
of one and two channels. In section V we extend the 
findings to the case of many channels. In section VI we 
come back to the A (3872) and comment on its decay to 
J/tp plus two and three pions, in section VII we comment 
on the independence of the results with the choice of the 
potential and in section VIII we outline our conclusions. 



II. THE X (3872) WITHIN COUPLED 
CHANNELS D°D*° AND D + D*~ 

In the A (3872) was plausibly explained as an I = 
dynamically generated state in coupled channels with 
positive C-parity. However, in that work the charged 
and neutral D mesons were put with the same mass, in 
which case one had a good isospin symmetry. In [l3| the 
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masses were taken different and a small isospin break- 
ing was produced. Summarizing the approach of [l2[ we 
call channels 1 and 2 the D°D*° and D+D*-. It was 
found in this work, using the hidden gauge Lagrangians 
adapted to the ST/ (4) flavor symmetry, containing ex- 
plicit breaking of the symmetry, that the potential in 
coupled channels, in s-wave, was very close to the type 



V 



FT 



FT FT 

V V 

FT FT 

V V 



(2) 



the label FT standing for field theoretical approach, and 
v FT is, in principle, a function of the invariant mass s 
(see Eqs. (7) and (12) of Ref. QJ]). To describe the dy- 
namics of the V(3872), which is placed quite close to the 
DD* threshold, it is sufficient to take the potential given 
in Eq. (7) of Ref. [l3j]) at threshold, neglecting in the po- 
tential all isospin breaking corrections induced by the dif- 
ference of masses between charged and neutral mesonifl 



FT 



m D m D , 



(3) 



with fjj ~ 165 MeV, the D— meson decay constant, and 
to jj and ttid* averages of the neutral and charged D and 
D* meson masses, respectively. The above interaction, 
should be considered with an ultraviolet cutoff in mo- 
mentum space of natural size for hadron interactions, 
A < 1 GcV. In numerical calculations we use mpo = 1865 
MeV, m B » = 2007 MeV, m D + = m D - = 1870 MeV and 
m D ,+ = m D »- = 2010 MeV. 

The Bethe-Salpeter equation in coupled channels in the 
on-shell factorization approach stemming from the use of 
the Nf D method [30-[32[ is given by: 

T FT - (1 - F FT G FT )"V FT (4) 

where both V FT and T FT are on-shel0, and G FT is the 
diagonal loop function for the two intermediate D and 
D* meson propagators. In the particular case of Eq. ((2J, 
Eq. (dJ) is trivially written as 



T 



FT 



V 



FT 



1 - uFTGfT - V ^G. 



FT 

22 



(6) 



We will assume that the T FT -matrix develops a pole, 
in the first Riemann sheet and below both thresholds, 
for the V(3872) resonance and the couplings g FT to the 



1 We will keep those isospin breaking corrections in the loop func- 
tion G that will be introduced below. As it was discussed at 
length in Ref. [T3l . they turn out to be quite relevant. 

2 The normalization is fixed thanks to the relation between the 
scattering matrix and the differential center of mass cross section, 
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FIG. 1. X decay mechanism into J/ip^(p) assuming the X 
to be a DD* molecule. 



channels DD* are defined such that in the vicinity of the 



pole 
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(7) 



with sr the squared mass of the resonance, which allows 
one to obtain the couplings via 



FT FT r / 

9i 9j = lim (s 

FTi2 



sr)Z 



FT 



-1 



f s (Gf? + Gf 2 T ) 



= [p ] 



(8) 
(9) 



All the couplings are equal in this case and, as we see in 
Eq. ((9|), they are independent of v FT and only depend on 
the derivative of the loop function G FT . One is tempted 
to interpret the couplings as the components of the wave 
function, and since the couplings are equal we would have 



o n*0\ 



|X(3872)) cx \D°D 



\D + D* 



(10) 



which represents a pure isospin 1 = state. Such an 
interpretation has some basis since the equality of the 
couplings is what makes the state behave as an I = 
state in the field theoretical approach. Indeed, think of 
the ratio of Eq. ([T]). The two pion and the three pion 
states in the decay of the A" (3872) correspond to a p and 
an uj respectively, according to experiment [33l l3~Ij . In a 
field theoretical approach the mechanism for this decay 
is depicted in Figjl] and the ratio for p and u decays 
would be given bjo 



Rnli,, — 



„FT,^fFT _FT/"iFT 
9l Wl ~ 92 U 2 2 
^FT^FT , n FT/~tf 
9l Wl + 92 U 2 



(11) 



The plus and minus signs in the numerator and denom- 
inator of Eq. (ITT]) are simply a consequence of the fact 
that J/ipp has 1=1 while J/ipuj has 1 = 0. If we had 
equal masses for the charged and neutral D mesons, the 
numerator of Eq. pT|) vanishes and the decay X — > J/ipp 
is forbidden, since it violates isospin. If the masses of the 
two channels are different, even taking the two couplings 



3 Technically, there can be a multiplicative coefficient in Eq. JTT}, 
but in the model used in ll.'l it is equal to 1. 
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FIG. 2. Loop function G FT for D° D*° and D + D*~ from 

El. 



A. The Lippmann Schwinger equation 

We need a potential V and to illustrate our results, 
which are general, we can take some easy form for it 
(other forms will be analyzed in section |VIT|) . We choose 
a separable function in momentum space with the mod- 
ulating factor being a simple step function, O. Thus our 
potential, already projected in s— wave, is assumed to 

(p'\V\p) =V(p',p) =vQ(A- p)0(A - p') (12) 

where p and p' stand for \p\ and \p'\ and A is a cutoff in 
momentum space. 

Let the Hamiltonian be H — Hq + V with Ho the free 
Hamiltonian. In this case the non-relativistic Lippmann 
Schwinger equation can be written as 



T = V + V 



E-Hn 



-T 



(13) 



or also as 



equal, as we found in Eq. ©, the numerator of Eq. (fTTj) 
does not vanish due to the difference in the loop function 
for each channel, see Fig. [21 and the decay X —> J/ipp 
is allowed. One could interpret this by saying that the 
X(3872) is an I — state but the intermediate loops 
in the decay violate isospin. However, the free Hamilto- 
nian (including masses) of the D°D*° and D + D*~ sys- 
tem does not commute with isospin, which implies that 
the X(3872) does not have a well defined isospin. The 
decays of the X(3872) resonance can provide informa- 
tion on this isospin mixture. In this work we will look 
at its J/ipp and J/ipuj decays, assuming transition oper- 
ators of zero (short) range. Within this scheme, we will 
show that the isospin violation in these decays is linked to 
the different probability amplitudes of finding the D°D*° 
or D + D*~ meson pairs, which form the X molecule, at 
short relative distances. This is the same as stating that 
the D°D*° and D + D*~ wave functions around the ori- 
gin will determine the isospin violations for these decays. 
Hence the couplings should not be understood as a mea- 
sure of the wave function components. We shall come 
back to this issue in what follows. 



III. COUPLING AND WAVE FUNCTION IN 
THE ONE CHANNEL CASE 



T = V + V- 



1 



E-H 

Taking Eq. (fTB"]) we can write: 
(P\T\P') = (P\V\P') + 

& 

fc<A E — mi — mi 



V 



(14) 



r7(k\T\p') (15) 



_ 

2/i 



where p is the reduced mass of the two particles that 
interact [1/p — \jm\ + 1/7712]. 
Eq. (fl~5j) has solution 



(p\T\p') = e(A-p)e(A- P ')t (i6) 

which can also be seen from Eq. (fT4"]l . with t given by 



t = v + vGt, 
G = [ d 3 P 



t 



1-vG 



1 



E - mi - m 2 - f^j 



(17) 
(18) 



We can see that Eq. (fl7|) is like the on-shell factorized 
equation (no integral left) of Eq. (gj), and G is indeed, 
up to a factoid, the non-relativistic reduction of the loop 



We will first study the non-relativistic dynamics of a 
bound state generated by the interaction of two particles 
of masses mi and 11J2, respectively. 



4 We use normalization (p\x) = e , P' x /(2tt)2 which means 
J d 3 x\x)(x\ = J d 3 p|p)(p| = 1 so that for a local potential we 

have <p"|V|p> = J d 3 xe- i< -P'-P^ s V(x)/(2Tr) 3 . 

5 The non-relativistic reduction of the scattering matrix, potential 
and two particle propagator loop function introduced in previous 
Sect. [IT] are related to those denned here by 



: 32lT 3 fJ,y/s V, 



: 327T 3 ^ V / S t 



(19) 
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function of two particle propagators regularized with a 
cutoff A, as is also usually done in the studies of hadron 
interactions with the on-shell Bethe Salpeter equation 
[351 . The name Bethe Salpeter equation was adopted 
in [35[ because there, relativistic meson propagators are 
used. For very weak bindings, it is sufficient to use the 
Lippmann Schwinger equation, which is what we do here. 
The poles of Eq. (|T7| occur for 



1 - vG = 



(22) 



which will occur for some value E a < mi + m 2 where we 
have a bound state. If the energy of the state is known, 
the above equation fixes the cut-off A, or conversely, if 
the cut-off is fixed, one can predict the energy of the 
state. 

In the work of [l2| the theoretical approach is based 
on the underlying hidden gauge formalism 36-38]. One 
assumes that it provides the potential V within a scale 
of momenta which is of the order of the cut off assumed 
here. Yet, one can assume certain uncertainties in V 
which can be compensated with changes in the range A 
in order to obtain the proper binding, as we show below. 
Indeed, using Eq. (|22|) . invariance under renormalization 
leads to 



where Eq. (fT7|) has been used for t and the l'Hopital's 
rule has been applied in the second equation. 

The integral for the G function defined in Eq. (fTBf can 
be performed analytically and we obtain 



(i{E n ) = - N//7T ( A — 7arctan ^~ 
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(27) 
(28) 



where < Eg =mj + m 2 — E a is the binding energy of 
the state a. The above equation allows us, using Eq. ([2"6"1) . 
to write the coupling g as 



9 2 = 



7 



87r/i 2 ( arctan 



7A 

7 2 +A 2 



(29) 



We note that g has a very smooth dependence on the 
cutoff A, which can even be removed (A — > 00). The 
coupling is mostly determined by the binding energy (see 
Eq. (fSTj) below) , specially in the limit in which the latter 
one is much smaller than the cutoff, and hence g is, to 
great extent, renormalization scheme independent. 



The wave function 



dK^ = ° 



(23) 
(24) 



1 dv __ 1 dG 
vdA ~ ~GdA 

and hence for a fixed E a , the bare potential should have 
a dependence on the cutoff to compensate that of the 
loop function G. 

B. The couplings of the state 

The coupling in this case is defined as g such that in 
the vicinity of the pole the scattering matrix behaves as 



t = 



E — E a 



and hence 



lim (E-E a )t 



G 



dG 
dE 



E=E„ 



327r 3 /iy / s 



with 



d 4 q 



(2tt) 4 q 2 -mf+ it (P - q) 2 - mi. + \ 



(25) 



(26) 



(20) 



(21) 



The Schrodingcr equation is given by: 

H\i/>) =E\4>) 



(30) 



where ip is an eigenfunction of H, the full Hamiltonian. 
We can write: 



(H + V)\1)) = E\1>) 

W) = 1 



E-H 



V\1>) 



(31) 
(32) 



which has the solution 

(p|V) = fd 3 k [ d 3 k'(p\ 



1 



E — Hq 



\k') 



x (k'\V\k)(k\t(>) 

e(A-p) 



E — n%i — n%2 — Jk<h 



(33) 

d 3 k(k\tp) (34) 



which gives us the wave function. Integrating Eq. (|34l 
over d 3 p, we obtain 



(35) 



1 - vG{E) = 



which is the condition to find the pole given in Eq. 

In Eq. (1341) , we determined the state wave function up 
to a constant, J fe<A d 3 k(k\^), which can be fixed from the 
normalization condition. Let E a < mi + m 2 be the solu- 
tion of the above quantization equation, its wave function 
will satisfy 



and P U P^ = s. 



d 3 p\(p\ 



1 



(36) 
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Note, that the wave function can be normalized because 
we are dealing with a bound state whose energy is below 
mi + TO2. From the above equation, one easily finds 



1 = ^ 



d 3 p 



1 



p<A 



E a - mi - m-2 - §^ y 



d 3 k(k\ip) 



>k<A 

and hence, it follows 



i! / d d k(k\ip) 

k<A 



dG 
dE 



(37) 



(38) 



E=E n 



We can now use the form of Eq. (fT4")) to solve the T 
matrix. We would have 



T = V+ V\m)(r 



1 



'E-H 



\m')(m'\V (39) 



where |m) and \m!) are complete sets of eigenstates of 
H . In the vicinity of the pole at E = E a we care only 
for the contribution of channel a, 



(p\T\p') ~ (i?|V 



1 



£7 — £? a 
= I d 3 kj d 3 k'(p\V\k)(k\a) 

' (a|fc')(*'W) 



E — E a 
v Ik<A d 3 k(k\a) 



-9(A-p)6(A-p') (40) 



E — E a 

exhibiting the form of Eq. (1161) from where we defined t. 
We can obtain the residue of t as 



g 2 = Urn (E — E a )T 

hi — Yhta. 



:■> / dk{k\a) 

k<A 



dG 
dE 



(41) 



E=E n 



where we have used Eq. (|38|) to get the last equality (a 
stands for ip here) and we note that this is the same result 
as in Eq. p6]l . 

The wave function in coordinate space can be equally 
evaluated: 



x\ip) = / d 3 p{x\p){p\4>) 



d 3 p 



(2tt) 3 / 2 



(42) 



Using Eq. §Q we fincfl 



(x\ip) = g 



2 1. 



Im / dpp- 



= ipr 



- mi - m 2 



£1 
2/j 



For large values of r this function goes as 

A 



4irr 



where 



A= -2(igV2~Tr{l + 0(l/A)} 



(43) 

(44) 

(45) 
(46) 



The exponential fall-off at large distances is controlled 
by the binding energy of the state a and the coupling 
of the state. This behavior follows the general rule of 
bound states outside the interaction region and, as we 
see, is largely independent of the cut-off A. For the same 
reason, A carries information on the interaction region. 



D. The meaning of the coupling in terms of wave 
functions 



By means of Eq. (|4"Tj) we know that (assuming a real 
wave-function for the bound state) 



d A k( 



k<A 



(47) 



and from Eq. (|4"2"|) we can obtain the value of the wave 
function at the origin in coordinate space: 



(x = 6\ip) = ip0) 



(2tt) 3 /' 



■m) (48) 



and so 



(2n) 3 ^G- 1 (EM6) 



(49) 



where we have also used the condition for the bound state 
1 — vG = 0. Since g hardly depends on the cutoff A, the 
wave function at the origin inherits the linear dependence 
on A exhibited by G(E a ) in Eq. (|2T|) . Yet, if v is known, 
the binding energy also fixes G(E a ) from the condition 
vG(E a ) = 1. 
Now we define 



rfj = gG(E a ) = {2^f' 2 ^{Q) 



(50) 



This constant will appear often in what follows. This 
is an important result concerning the problem at stake: 
the coupling, up to the factor G(E a ) 7 is a measure of the 
wave function in coordinate space at the origin. 



It is possible to write Eq. (143 I t in terms of the analytical functions 
sine integral and cosine integral. 
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E. The limit of small bindings 



IV. TWO COUPLED CHANNELS 



Taking the limit for small values of 7 in Eq. (|29|) we 
see that 



r 2 7 
hm g = - „ „ 
7^0 * An 2 fi 2 



(51) 



a result well known (up to a normalization depending on 
definitions) [HS|]. 

This result can also be obtained from the general form 
of the scattering amplitude at low energies. Indeed, let 
us recall the form of the for s-wave scattering amplitude, 
/, close but above threshold, 



/-!(£) = fccoU- ifc 



1 r k 2 

h 

a 2 



-ifc, (52) 



with d the phase shifts, a and tq effective range param- 
eters and fc = y/2fi(E — mi — m®). The analytic contin- 
uation below threshold to the energy of the bound state 
E = E a reads (we assume E a is close to threshold) 



roT 



+ ••■7 



(53) 



where we have taken k(E a ) = i"f. The inverse of the scat- 
tering matrix must vanish at E = E a , as it corresponds 
to a bound state, and the limit: 



jjj 

f = lim (E - E a )f(E) = lim - „ ° 
-1 7 



2 



7 



0(Y) (54) 



E—E a 



which agrees with Eq. (50) of Ref. |42i |. 

This result is equivalent to the one in Eq. (I5T1) up to a 
normalization which is easy to get recalling that 



f 



8tt, 



-Air 2 [it 



(55) 



As a consequence the coupling that we are using becomes 



9 



7 



4tt 2 



-0( 7 2 



(56) 



which is the result obtained in Eq. (f5Tj) . 

A final remark concerns the comparison of this result 
with the couplings defined in 13[ and in general in studies 
using the chiral unitary approach [35[ where G is defined 
in a field theoretical approach in terms of two relativistic 
propagators (Eq. ([2"T])). Note that from Eqs. (TIT)]) and 

(da 



FTFT 



G tl v 



Gv 



(57) 



which guaranties that the position of the pole remains 
unchanged, since it is determined by the condition Gv = 
1. Besides, from Eq. (fT9|). we trivially find 



FT 



(647r 3 /0 3 g 

E a (16tt 7 /tf ( 7 ->0) 



(58) 
(59) 



A. The couplings 

We work out in this section the two channel problem 
for the particular case of the A(3872) using a dynam- 
ics determined by the potential of Eq. ((2J, which is an 
isoscalar operator (i.e., it is diagonal in the isospin ba- 
sis). We will work first within the Quantum Mechanics 
formalism that is adequate here, since the mass of the 
A (3872) resonance and that of the charged and neutral 
DD* pairs differs in just few MeV. We will use again a 
cut-off A in momentum space, and the 2x2 matrices T 
and V will encode step functions 



(p'\V\p) =V(p',p) 
(p\T\p') =T(p',p) 



u6(A 

te(A- 



-p)0(A 
p)0(A- 



-P) 

V) 



with 



v v 



(60) 



(61) 



The Lippmann Schwinger equation in the coupled chan- 
nel space reads 



t = (1 - vG)-^ 
1 



1 - vGn - vG 



22 



where 



G = 



G11 
Go, 



Ga — 



d 3 p 



p <a E-Mi-% 



(62) 
(63) 

(64) 



2/i, 



with E the relative energy including the mass of the par- 
ticles and Mi and M2 the thresholds of each channel 

Mi = m D a + rrifj.o, Ah = m D + +m D ,- (65) 

and /ii and \ii the reduced masses of the D°D*° 
and D + D*~ systems respectiveljQ- Assuming that the 



7 The correspondence of the Quantum Field Theory and Quantum 
Mechanics, in the nonrelativistic limit, scattering matrix, two 
particle propagator and potential matrices reads 



-iFT 



1 



fi 2 G n 2 



where 



327r 3 yi 



(66) 

(67) 
(68) 



(69) 



and ^t - 2 the inverse of the above matrix. The Bethe Salpeter 
equation implies, in the non-relativistic limit, the Lippmann 
Schwinger equation J62}. Thus, it trivially follows 



det(l - V FT G FT ) = det(l — v G) 



(70) 
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X(3872) is a bound state of the system, its mass (E a < 
Mi < M 2 ) will be obtained by requiring that the denomi- 
nator of Eq. (|63|) will vanish (pole of the t— matrix in the 
first Riemann sheet and below all thresholds). Clearly 
when E a — > Ml the D°D*° is loosely bound and the 
D + D*~ is bound by about 8 MeV. The explicit expres- 
sions for Gu(E a ) are given by Eq. (|27p with 



7» 



E a . 



(72) 
(73) 



Let us now pay attention to the couplings of the bound 
state to each of the two channels. Since all elements of 
the matrix v are equal, both couplings g\ and g 2 are the 
same: 



9l 



lim (E — E a )Uj 
dG 22 



E—tE, 

(dG n 



V dE 



dE 



E=E n 



On the other hand, as in Eq. (|58|) , we obtain 



FT 



!1 



= (64n 3 flEl) 



fJ 



(74) 



(75) 



with /i, the average of the fj,\ and [i 2 reduced masses, as 
discussed in the footnote [7] By using Eqs. (|74|) and (JT5J) , 
we find not only a qualitative, but also a quantitative 
good agreement with the results for the coupling shown 
in Fig. 3 of Ref. [l3[ . For instance, if the neutral channel 
is bound by 1 (0.1) MeV, the denominator of Eq. ([55)1 
will vanish for a value of the cutoff A of around 680 (653) 
MeV, with v given by Eq. (1711) . This leads to a cou- 
pling g FT of the order of 5400 (3200) MeV, in reasonable 
agreement with the result in Ref. |13| . 

In the limit when E B1 — > 0, we have ^j E x \ E _ E — > 00 
and we find 



9l 



7i 



47T 2 /if 



(76) 



thus, both couplings go to zero as ^/tT 



which guaranties that poles are placed in the same position in 
both approaches. Finally, note that a potential in the field the- 
ory approach of the form in Eq. J5J does not lead to a quantum 
mechanics potential of the form assumed in Eq. 1611 . unless that 
both reduced masses fii and fj.2 are taken to be equal in Eq. ((69} . 
This is an excellent approximation for the case of the coupled 
channels D°D*° and D+D*~, and we have done so here when 
relating Quantum Field Theory and Quantum Mechanics quan- 
tities. Thus, we have taken /ii ~ fi2 ~ p., with p, some average 
reduce mass. Hence, the Quantum Mechanics potential, v, near 
threshold and for an ultraviolet cutoff of natural size for hadron 
interactions A < 1 GeV, can be now approximated by 



1 



327T3/2 



(71) 



as deduced from Eq. (J3J , with the further approximation fJ.y/s ; 
m D m D * . 



B. The wave function 



For the bound state we have now a two component 
wave function, representing each of the D°D*° and 
£)+£)*- channels, and the Schrodinger equation reads 



(H + VM=E\iP) 



The solution to this equation is given by 
1 








(77) 
(78) 

(79) 



(P IV') = t 

E-M 2 -p' 2 /2n 2 

x / d 3 k(p\V\k)(k\ip) (80) 



which represents two coupled channels equations 
Q(A-p) 



(p\ipi) = v 



E — Mi 



(p \tp 2 ) = v 



x / d :i k ( (fc|Vi) + (fc|V 2 ) 

k<A K 

e(A- P ) 



(81) 



x / d\ Uk^x) + (k\ip 2 )) (82) 

which require to know J fc<A d 3 k(k\^i) for its solution. To 
evaluate this latter magnitude let us integrate over p in 
Eq. (JHSJ and we get 

d 3 p(p\i>) =Gv [ d 3 p(p\^) (83) 

p<A Jp<A 

an algebraic equation that requires for its solution 

det(l — Gv) =1 — vG u - vG 22 = (84) 

This equation is satisfied for the poles, E — E a , of the 
t matrix corresponding to bound states (see Eq. ([6"3"|)1. 
Eq. ([55)1 can be now be solved and we obtain 



p<A 



1 — V U 22 J p< j\ 



u ll Jp<A 



(85) 



where G% = G U (E = E a ). Thus, Eqs. flU) and ^S2J can 
be written as 



(pIV'i 



i e(A - p) 

GfiE a ~M 1 ~£ I j, \ 



d 3 fc(A#i) (86) 



(PW = ^ - P) p - 2 / d 3 fcWi) (87) 

t*ii £ Q - M 2 - Jk<A 
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If we define the partial probability 

Pi = {^i) = J dM{p\^)\ 2 (88) 
we can further use the total normalization condition 

i = p 1 + p 2 ^ [ rfMl<pl^>l 2 + l(pWI 2 } 

Jp<A 



1 dG n 



1 dG 



22 



E=E a 



[Gfja dE 
d 3 /fc(fc|Vi) 

'fc<A 

from where, using Eq. (1741) 



[G^f dE 



E=E„ 



d 3 p{p\ipi 



p<A 



and hence 



Pi = -sr 



dG 



li 



dE 



E=E 



p ,,2., arctan 
r\ _ Mi72 



^2 = -g 2 

(A) - 7iA 



dG 2 



(89) 



(90) 



(91) 



E—E a 



f+A 2 



^2 M271 



^ [l + C(A _1 ) 
7i 



arctan 



(72) 



72 A 



^[^©(A- 1 )] (92) 



The neglected terms are finite range corrections, which 
in our case are represented by the finite cut-off. On the 
other hand, assuming real wave functions, Eq. to- 
gether with Eq. (|55T) lead to 



(93) 
(94) 



gG^= [ d 3 p{p\i> 

Jp<A 



gG 



22 



d 3 p(p\ip 2 



p<A 



The wave functions in coordinate space would be again 
given by means of Eq. (|43|) using fix for ipi and \X2 for ip2 
and substituting wi\ +ni2 by Mi for each component. At 
long distances this means, see Eq. (|44|) . 



(zlV-iL 



4irr 
A 2 



47rr 



(95) 
(96) 



where we have used Eq. (|73|) . Once again J p<A d 3 p(p\tpi) 

can be interpreted (up to a constant factor (2tt) i ) as the 
wave function at the origin, as done in Eq. (|4"5)l and then 
Eq. (|9"3")l and Eq. can be rewritten as 



gG^^(27r)^M0)=^i 
gG%2H^) 3/2 MQ)=ip2 



(97) 
(98) 




r[fm] 



FIG. 3. Neutral and charged Wave function components for 
a D°D*° binding energy of 0.1 MeV. In the upper panel one 
can see the value of the wave functions at the origin for both 
channels, in the lower panel we plot the probability density 
for each channel. 



^From the above expressions, we trivially find the ratio 
of wave functions at the origin 



4i 



Got 
22 

Ga 
11 



(vG a n y 



(99) 



and their departure from unity provides a tangible mea- 
sure of the isospin breaking in the interaction region. On 
the other hand, 



dA 



7T M2 72 ~ 7l 

2 /ii A 2 



+ C(1/A 4 ) (100) 



which shows that, though both ipi and tp2 depend greatly 
on the cutoff, such a dependence is much reduced in their 
ratio. 

In Fig. [3] we show the two wave function components 
for the problem solved in 13[ with only two channels 
and the binding energy for D°D*° of 0.1 MeV. The os- 
cillations are caused by the sharp cut-off A introduced to 
regularize the Lippmann-Schwinger equation. Eqs. (|97|) 
and (j9"8|) and the asymptotic behavior of the wave func- 
tions given by Eqs. (1951) and (|96[) clearly show that, in 



the limit of zero binding energy of the D°D*° compo- 
nent, this wave function extends up to infinity while the 
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D + D*~ component is restricted in space because of the 8 
MeV binding. The probability to have the D°D*° com- 
ponent becomes much larger than that of the D + D*~ 
component (see Eq. (j9"2")) ) and we could think of the 
A (3872) as a D°D*° molecule. While technically correct 
from the point of view of probabilities, this interpretation 
is misleading concerning physical processes, like decays, 
because these require Hamiltonians of short range, zero 
range ordinarily in effective field theory, such that what 
matters in these processes is the wave function at the 
origin. For instance, what will determine the I — char- 
acter of the wave function will be the ipi(0) and ^(O) 
magnitudes, not the probability integrated in all coordi- 
nate space. Hence, Eqs. (|97|) and (|98l) indicate an isospin 
breaking, with respect to the 1 = combination, given 
by the differences between Gfi and G 22 (see Fig. [5]). 

Eq. (fTTj) . related to the decay process depicted in 
Fig. [TJ has now an intuitive representation to the light of 
Eqs. (HZ]) and (JM|)- The amplitude for 



is given by 



M 



X(3872) ->■ J/ipp{uj) 

g%GnFi + g 2 G 22 F 2 
(27r) 3 / 2 (V-i(0) J F 1 1 +^2(0)F 2 ) 



(101) 



(102) 



where F\ and F 2 are isospin factors for the vertices 
L>°L>*° -» J/ipp{u) and D+D*- J/tpp{uj), F 2 /F 1 = 
— 1 for the p and F 2 fF\ = 1 for the uj. Certainly this 
should be the case for many channels and we shall see 
the generalization in the next section. 

Eq. (I102p can alternatively be interpreted as 

M = J d 3 p ip\tpx)t D oO,o^j/^ p{ ^ 

+ J d 3 p(p\ip 2 )t D+D ,-^ J/ilp{u) (103) 

assuming that the range of the t amplitudes is very short 
compared to the extension of the wave functions, essen- 
tially that the t are constant functions in momentum 
space, as one has in field theory vertices stemming from 
a contact Lagrangian. 



A. The couplings 

We can write for the T matrix 



det(l - vG) 
where A is defined as 

A = [det{\ - vG)] {I -vGy 1 



(105) 



(106) 



This matrix is introduced to single out the source of the 
pole in all channels which is given by the condition 



det(l -vG)=0 



(107) 



We have now 



9i9j = lim (E - E a )Tij 



(Av) 



^det(l-vG) 
9j _ (Av) zj 



E=E n 



9i {Av)i 



(108) 
(109) 



E=E a 



We can see that gj / gi is a ratio of two matrix elements of 
matrices without singularities. This means that if gi — > 
as a consequence of having the binding in channel i go- 
ing to zero, then all the couplings to the other channels 
coupled to channel i will also go to zero. This is also 
obvious from Eq. (|108|) since g? — > because the denom- 
inator contains -j^- and, one has ^jff- — > 00 for all cases. 



B. Wave functions 

Eqs. (|8Tj) and (|82[) can be generalized as 

1 



(p|^) = e(A- P ) 



x y^vij I d a fc(fc|^-) 

lk<A 



E - M *-k 



(110) 



which upon integration leads to 



V. GENERALIZATION TO MANY CHANNELS 

We have now 

(p'\V\p) = v0(A — p)0(A — p') (104) 

where v is a N x N matrix with N the number of channels. 

The expressions that we have obtained can be gener- 
alized to many channels and we derive here some useful 
expressions. 



/ d 3 p(p |V0 = Gu ^2 v i. 
J P <A ,• 



d 3 k(k\ipj 



fc<A 



which in the language of Eq. (|9T)l and Eq. 

ipi = Gu y^Vjj^j 

j 

tp = Gvip 



(111) 
reads 

(112) 



where Eq. (|112p is written in matrix form and it requires 
that det(l — vG) = for its solution, which is guaranteed 
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for the bound eigenstate. Eq. (|112l) can also be written 
as (G- 1 =G- 1 (E a )) 



G- 1 ^ = vi> 



(113) 



which allows to rewrite the equation for the wave func- 
tions Eq. (jllOl) as 



diag 



O(A-p) 



E,, 



M, 



(H4) 



which gives the wave function in momentum space in 
terms of the wave function in coordinate space at the 
origin. These equations are the generalization of Eq. (|86j) 
and Eq. {57]) together with Eq. (|55|). 

Let us now use the normalization condition 



i 



x - dG... 



E=E a 



(115) 



We can now take advantage of Eq. (|114[) to define the 
couplings in terms of the i/ji. For this we use the ver- 
sion of Eq. (|14p for the Lippmann Schwinger equation, 
recalling that close to the pole of the eigenfunction of 
the Hamiltonian, associated to the energy E ai only 
this state contributes in the sum over eigenstates of i?, 
and we find 



Tij — Vij 



E 



i 



/ d k(k\ip m ) 

k<A 



E-E n 



d 3 k'(k'\^ n )v ri 



(116) 



a Jk'<A 

which means that 

9i9j = Vimi'mVnji'n 



E=E a 



(117) 



from where we conclude that 



9iG% = & 



(118) 



as we found in Eq. (p?T)) and Eq. (f9"5)) in the two channel 
problem. This allows to reinterpret Eq. (I115P in terms of 
the couplings and we find 



dG, 



\ 2 



dE 



■1 



(119) 



E—E a 



which is the generalization of Eq. (f74"f . 

Eq. (| 1 19[) is interesting because when one channel be- 
comes loosely bound then the loop derivative for this 
channel goes to infinity while the other derivatives re- 
main finite. In this limit wc get, if channel 1 is loosely 
bound 



lim of 



2 dG u 



dE 



1 



E—E a 



dGu 
dE 



(120) 



E=E a 



which is the same result obtained for one channel in 
Eq. (jUJ). Thus, in this limit the coupling of the loosely 
bound state goes to zero as the binding energy goes to 
zero. On the other hand, Eq. (|109j) guarantees that all 
the other couplings will also go to zero since the matrix 
Av is not singular. This result was also found in [42[ 
although derived in a different way. 



VI. DECAY WIDTH OF THE X(3872) 

After these clarifications we would like to go back to 
the ratio of the decay width of the A(3872) to J/tpp and 
J/i/jlj of Eq. (TTTj) . As discussed in [l]| , the ratio of widths 
was given by the square of Eq. (fTT|) times the factor to 
correct for the phase space of p decaying to two pions 
and the w decaying to three pions: 



Ga 
11 



U 22 



B{X -)• J/^tttttt) V G n + G 22 



J^qS (s,m p ,T p ) 9 (m x 



yfs) ds Bp 



J °° qS (s, rn w , r w ) 6 (m x - mj/^ - y/s) ds 



(121) 



where B p and B p are the branching fractions of p decay- 
ing into two pions (~ 100 %) and u> decaying into three 
pions (~ 89 %), q is the center of mass momentum of the 
outgoing meson pair in each channel and value of s, and 
S (s,m,T) is the spectral function of the mesons given 
by: 



>S (s, m, r) 



7T 



1 



iTm 



(122) 



In [l3j it was found, using dimensional regularization for 
the loops, 



B(X -> J/i>ir + Tr-Tr ] 

B(X -> J/^TT+TT") 



1.4 



(123) 



which is compatible with the experimental value 1.0 ±0.4 
from 0. 

Now let us assume that we take seriously that there 
is only one channel, the D°D*°. Then the ratio of the 
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widths is 



R 



(D°D'°) _ ( % l ) l t D Q D- c >^J/^f, 



p/u 



lplt D o 



= 1 



(124) 



which is about 30 times bigger than the value obtained 
for this ratio (0.032) in (l3lp[. When we take into account 
the phase space for the decay into p and u> and the p and 
lu branching ratios into two and three pions, with the 
ratio in Eq. (|124l) we find 



B(X -> J/lpTT+TT-TT ) 



= 0.05 



(126) 



which is about a factor 20 times smaller than experiment. 

It is thus clear that the charged components of the 
wave function have played an essential role bringing this 
branching ratio close to experiment and this stresses that 
the wave functions at the origin for each channel, and not 
the probabilities of finding the state in a single channel 
alone, is what determines the isospin nature of the state 
in coupled channels. Indeed, the X(3872) wave function 
would read 



:r\r = ^ Xl=0 



V2 
^/2 



Xi=i 



(127) 



with X/=o,i scalar and vector isospin wave function 
spinors. In the charge basis used in Section IV, the 
isospin wave functions are 



Xi=o 



1 

71 



Xi=i = 



1 

71 



i 

-i 



(128) 



As mentioned in the introduction, the A (3872) does 
not have well defined isospin because the free Hamil- 
tonian (including masses) of the D°D*° and D + D*~ 
system does not commute with isospin. The mixing 
depends on the relative distance r between the pseu- 
doscalar and vector mesons. Thus with transition op- 
erators of zero range, one easily understands that the 
ratio of branching fractions of Eq. (JT|) is determined by 



the ratio 



(i - vWv>i) / (i + vv^i) 

gives the ratio of isospin 1 to isospin probabilities at 
r = 0. As we shall see in the next section (see R p / U of 
Table H]), this ratio is of the order of 2%. 



which in turn 



8 Note that 



( Gil 


~ G 22 \ 


' l — 1p2 /lf>l \ 






I 1 + ^2/^1 / 



(125) 



and thanks to Eq. I|100| l the ratio of wave functions at the origin 
depends little on the ultraviolet cutoff A. 



VII. RESULTS FOR OTHER FORMS OF THE 
POTENTIAL 

One might think that the results obtained are specific 
of the type of the potential chosen in Eq. (jSCJ)) . but the 
results are actually very general. To show that this is the 
case, we use other potentials. Let us consider a separable 
potential where we substitute the sharp cut off by a form 
factor 



(p'\V\p) = vf(p)f(p') 



(129) 



where v is a N x N matrix with N the number of channels. 
The results of Sect. V follow nearly identically substitut- 
ing the 6(A -p) by f(p). Eq. (fT05>Eq. ([109} are the 
same, but G is now given by 



d d P f 2 (p) 



1 



(130) 



and the wave functions are now given by 

1 



(m) = f(p) 



E a -Mi-$- 



X / d 3 kf(k)(k\^ 



Eq. CLID-Eq. (fTT51) follow, but now 
which allows to write Eq. (|131[) as 

m 



(131) 



(132) 



(jp\ip) = diag 



(133) 



and again we find Eq. (|118l) 



.9, - WGg 



(134) 



and Eq. (ITT^)) -Eq. dT2U| also follow. 

Everything is identical as before, but now ip is not, up 
to a factor (2ir)2 , the wave function at the origin (it 
would be if we removed f(p) from Eq. (|132[1 ). To see 
the meaning of ip we write f(p) in terms of its Fourier 
Transform 



1 



(2tt) 3 / 2 



d 3 xf(x)e l P s 



(135) 



and the wave function of Eq. (I133P also in terms of its 
Fourier Transform 



1 



(2tt) 3 / 2 



d 3 xe-^ s iPi(x) (136) 



Then upon integrating explicitly over k in Eq. (|132l) we 
find 



(137) 



ipi = / d 3 x^(x)f(x) 
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We performed explicit calculations using a gaussian form 
for f(p) 



fip) = e 5 



3 „-■ iafA' 



(138) 



and a Lorentz form 



m 



A 2 



A 2 +p 2 



\x\A 



2 \x\ 



(139) 



We can see that f(x) has a range of 1/A ~ 0.2 — 0.3 
fm, a range much smaller than the extension of the wave 
function. Thus ipi gives the average of the wave function 
in the vicinity of the origin, while in the case of the sharp 
cut off one finds exactly the wave function at the origin. 

One sees again that this average value of the wave func- 
tion at the origin is what governs the decay process of the 
resonance. Indeed, if we go to Eq. (|103l) to get the am- 
plitudes for the decay mechanism of Fig. [T] we would 
find: 



M= J dp (p\tpi)f(p) t D <>D*°^J/4p{u) 

+ / d 3 p(p\iP 2 )f(p)t D+D «-^ J/i!piuj} (140) 



where we have taken for the DD* —> 
amplitude, for consistency, f(p)t D jj 
constant. Eq. (|140[) reads 



J/ipp(u>) transition 
--►J/iMw) with 1 a 



as we obtained before, and according to Eq. f|134[) now 
9iGa — ipi. 

It is interesting to compare the results with the form 
factor with those we had before. We proceed as follows, 
we take the same strength for the potential as before and 
determine A to get the neutral channel bound D°D*° at 
the same energy. In Table Q] we show the values of g, 
ip and A for the three approaches: the sharp cut off, a 
gaussian form factor and a lorentzian form factor. We 
see that the differences for g and ipi, the relevant mag- 
nitudes in the X (3872) decay, are small. Note that the 
value of the R p / U of [l3[ was obtained using dimensional 
regularization and zero binding energy for the D°D*° 
channel. If the calculation is made in the same scheme 
but with Eb— 0.1 MeV one obtains i? p / w =0.025, which 
differs from the values on Table U The bulk of this dif- 
ference can be attributed to differences between using 
hard cut off in momentum space or dimensional regular- 
ization. Indeed, using the relativistic treatment of [l3j . 
but using a hard cut off, provides a value for i? p / w =0.020 
at E B =0.1 MeV. 

We also plot in Fig. @]the wave functions obtained with 
the three approaches and, as we see, they hardly differ 



above 0.5 — 0.6 fm. The differences at the origin cancel in 
ipi, where a smearing around the origin is involved, and 
particularly in the ratios tpi/ip2- 

The robustness of the results against changes in the 
form factor suggests that a large cut-off approximation 
should work fairly well. Indeed, in this limit all cut- 
off procedures merge into a single one, and taking into 
account that Gu{E a ) — Gu(Mi) is finite we get 



Gu(E a 



Sirpi 



+ 0(A- 1 ) 



dp[f(p)} 2 ~ hi 
o z 



(142) 



which shows that what matters is the integrated strength 
of \f(p)] 2 and corresponds to using a common subtraction 
constant with a different cut-off interpretation, 



dp[f(p)f 



A 



Sharp 



-A 



Gauss 



-A 



Lorentz 



(143) 



This identification provides a simple rule relating the dif- 
ferent cut-offs Asharp, A Ga uss and A Lorcntz which works 
very well as can be checked from Table HI Keeping the 
leading terms in Eq. (|142[) and taking \i\ = fj,% = p, we 
obtain the following remarkably simple analytical results 
by using the bound state condition, Eq. (|84l) . 



(144) 
(145) 



' dp[f{p)?=*-{ ll+l2 )+ 2 * 2pD 



R 



p/u 



(7i ~ 72) V 

64/47T 2 



(141) which yields 



dp[f{p)} 2 = 665MeV, R p/u = 0.025 (146) 



Next we want to connect the present results with those 
obtained in [T3, EH in the relativistic approach. Actually, 
by matching the relativistic one-loop integral calculated 
within dimensional regularization with scale v (called fi 
in [12| | ) with the non-relativistic propagator in the heavy 
meson limit tod,wd* 3> 71,72 (see Eqs. ([20]) . ([2Tj) and 
(US])), we ge10 



dp[f(p)f 



(mo + TTlD* )an 



+ m D log(-£) + m D , log(^fl)(147) 



where an is a dimensionless subtraction constant which 
depends on the scale v. For v = 1.5 GeV the subtraction 



9 We are using averaged values of masses for D and D* in this 
formula. This is appropriate since no differences of masses appear 
in the formula. This allows the matching between relativistic and 
non-relativistic expressions. 
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TABLE I. Comparative results for different potentials for a 
D°D*° binding energy of 0.1 MeV. 



Form 
Factor 


A 

[MeV] 


g FT 
[MeV] 


V;i(o)/V 2 (o) 




V>i 


4>2 




Sharp 


653 


3202 


1.31 


1.31 


3.29 


2.50 


0.018 


Gauss 


731 


3238 


1.20 


1.29 


3.30 


2.56 


0.016 


Lorentz 


834 


3254 


1.17 


1.28 


3.31 


2.58 


0.015 



Sharp Cut off 




1 2 3 4 5 6 
r[fm] 

FIG. 4. Wave functions for different form factors in the po- 
tential. 



constant used in Ref. [13] was an = —1.185 yielding 
J °° dp[f(p)] 2 = 651 MeV in fairly good agreement with 
Table U (Asharp)- Neglecting the finite cut-off corrections 
as in Eq. (|142l) one obtains now R p / U = 0.026. 

VIII. CONCLUSIONS 

With a view to the structure of the X(3872) as a pos- 
sible coupled channel bound state of mostly the D°D*° 
and D + D*~ and other minor channels, we have stud- 
ied the meaning of the couplings, which one determines 
from the residues of the scattering matrix at the poles, 
in terms of wave functions for the different channels. We 
have done the study in one channel, then in two chan- 



nels suited to the AT (3872) resonance and then we have 
generalized the results to many channels. Interesting re- 
lationships are obtained which shed light on the field the- 
oretical approaches to reactions from the perspective of 
wave functions. Essentially we find that the couplings 
are proportional to the value of the wave function at the 
origin in coordinate space or the averaged value within 
the range of the interaction. They are not sensitive to 
the wave function at long distances which is governed by 
the binding energy, and we also find that this averaged 
value of the wave function at the origin, ip = gG, is the 
only information that is needed when dealing with short 
range processes, like those provided in terms of contact 
Lagrangians in field theory. We also found that the val- 
ues ipi were very stable against assumed shapes of the 
potential once the binding energy is fixed fulfilling the 
quantization condition det(l — VG) = 0. 

We also find that, when one channel becomes loosely 
bound, the couplings to all coupled channels go to zero. 
Even if in terms of probabilities the loosely bound chan- 
nel, whose wave function extends up to infinity, has the 
largest probability, what matters in the reactions is the 
averaged values of the wave functions at the origin that 
determine the dynamics of the processes and the under- 
lying symmetries like isospin. The isospin violation in 
particular is tied to the ratio of wave functions around 
the origin (for short range processes), which goes 

to a finite limit when the binding of the "01 component 
goes to zero. 

When coming to the AT (3872) case, which can corre- 
spond to the D°D*° channel very loosely bound and the 
£)+£)*- bound by about 8 MeV, we find that the wave 
functions at the origin for the two channels are simi- 
lar, suggesting that one has a state with 1 = 0, with 
small isospin breaking, even if the probability to find the 
D°D*° component in the full space is much larger than 
for the D + D*~ component. A precise measure of the 
isospin admixture is given by the ratio tpi/i/)?, which is 
very stable and has a value around 1.3, the value of 1 
corresponding to a pure 1 = state where the decay 
X —> J/4>p would be forbidden. 

The consideration of the charged D + D*~ component 
to describe physical processes is so important that if it 
is neglected one finds a ratio of ^x-^//^+ ■J'-Jo) twenty 
times bigger than experiment. 

The work done has also an academical component. 
Some useful expressions, as well as exact analytical so- 
lutions for the wave functions in coupled channels have 
been given. The work also shows a different perspective 
on the on-shell approach to the scattering matrix based 
on the N/D method used in all modern works of chiral 
dynamics in coupled channels, by means of which the cou- 
pled Bethe Salpeter integral equations become algebraic 
ones. The suitable choice of the potential in momentum 
space that we made gives rise to the same equations as 
in the field theoretical on-shell approach. The analytical 
expressions found can be very useful to give alternative 
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interpretations of results found in the chiral unitary ap- 
proach, or in general in unitary coupled channels meth- 
ods in many physical processes. 
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